Properties of the solutions of delocalised coagulation and inception 

problems with outflow boundaries 

Robert I. A. Patterson* 

February 19, 2015 


Abstract 


in 

t-H 

o 

(N 

<D 

O 

(N 


PH 

< 

a 


Well posedness is established for a family of equations modelling particle populations undergoing 
delocalised coagulation, advection, inflow and outflow in a externally specified velocity field. Very 
general particle types are allowed while the spatial domain is a bounded region of d-dimensional 
space for which every point lies on exactly one streamline associated with the velocity field. The 
problem is formulated as a semi-linear ODE in the Banach space of bounded measures on particle 
position and type space. A local Lipschitz property is established in total variation norm for the 
propagators (generalised semi-groups) associated with the problem and used to construct a Picard 
iteration that establishes local existence and global uniqueness for any initial condition. The unique 
weak solution is shown further to be a differentiable or at least bounded variation strong solution 
under smoothness assumptions on the parameters of the coagulation interaction. In the case of 
one spatial dimension strong differentiability is established even for coagulation parameters with 
a particular bounded variation structure in space. This one dimensional extension establishes the 
convergence of the simulation processes studied in [Patterson, Stock. Anal. Appl. 31, 2013] to a 
unique and differentiable limit. 

Keywords: Coagulation, advection, existence, uniqueness, regularity, Banach ODE, propagator, 
boundary 
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1 Introduction 

° i-i 

Smoluchowski 18| introduced equations for the concentrations of particles of different sizes undergoing 
ly^oagulation in a spatially homogeneous population. 

o 

<N 
O 
in 

where c(t, y ) is the concentration of particles of size y at time t and K is a symmetric function defining 
nthe ‘reaction’ rates. The Smoluchowski coagulation equations can be regarded as describing a system 
of binary reactions involving an infinite number of species, but with a very structured, although non- 
sparse set of rates and CD abstractly written c = R(c). The model therefore extends naturally to 
a reaction-transport problem for spatially inhomogeneous populations of coagulating particles of the 
general form c + Ac = R(c ) for some transport operator A. 

Since coagulation is a binary reaction in which every possible pair of particles may coagulate, 
the equations are, even in the spatially homogeneous case, non-linear and more significantly non¬ 
local in particle size (size may here be generalised to ‘type’). The first existence results for the 
Smoluchowski coagulation equation and its extensions were based on convergent sub-sequences of 
approximating stochastic processes. The first convergence result of this kind with simple diffusive 
transport of particles is due to Lang and Xanh {§], generalisations were achieved by Norris flfl. [Til] . 
Wells [19] and Yaghouti et al. [21]. This is quite a natural approach, because the equations are based on 
a microscopic stochastic model and related stochastic processes have also proved fruitful for numerical 
purposes going back to Marcus [91] and Gillespie [s]. 

The results just mentioned are essentially compactness results and say nothing about uniqueness 
of the limiting trajectories, much less of uniqueness for the solutions to the Smoluchowski equation 
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and its extensions. Convergence and uniqueness were proved together by Guia§ [6| who modelled 
diffusion as a random walk on a lattice and used a more functional analytic approach. Going further 
in this direction one is led to regard the Smoluchowski equation and its extensions as an ODE on 
a Banach space and to proceed via a locally Lipschitz source term and a Picard iteration method 
to show existence and uniqueness in some functional setting. The general strategy is presented in 
chapters 5&6 of 17]. Applications to Smoluchowski problems are given by ,20. 1,13:] and the works 
cited therein. 

Especially when approaching the Smoluchowski equation from the point of view of stochastic 
particle systems it is natural to think of measure valued solutions. A particle system is identified with 
its empirical measure and thus instead of functional solutions one is led to look at measure valued 
solutions in a weak setting. To give the concrete example that will be the focus of this work: A 
solution (with a given initial condition) is a flow of measures fit on positions in X and particle types 
(sizes and potentially additional details) in y satisfying 


d_ 

dt 


1 Xxy 


f(x,y)/M(dx,dy) 


/ u t {x) • V/(x, y)nt(dx, dy) + / f(x, y)I t {dx, dy) 

Ixxy Jxxy 

+ x I [f{xi,yi + V 2 ) - f(xi,yi) - f{x 2 , y 2 )\ 


'(Xxy) 2 


K(y 1 , y 2 )h(x 1 ,x 2 )m (daq, dyi)y t (dx 2 , d y 2 ) (2) 


for all / in a class of functions D to be specified below. Here the problem has been moved from the 
strong formulation of dTJ) to a weak setting; a transport operator ut • V (the dual of the A mentioned 
above) has been introduced and the delocalisation of the coagulation specified via a function h, which 
may be regarded as a modifier. A particle source term It has also been added, which is relevant for 
many real-world applications as discussed later. 

Signed measures can be regarded as Banach space under a wide range of norms and equation ([5]) 
interpreted as a Banach space valued ODE and Picard-like fixed point strategies introduced. An 
important insight of the monograph [ 7 ] was to exploit duality of linear operators and norms between 
measures and appropriate spaces of test functions in pursuit of this programme. In this way one 
performs most calculations for operators on test function spaces, which are a little easier to work with 
than operators on spaces of measures. Measure valued solutions are also the topic of 11], which also 
uses a linear operator approach, but uses approximation rather than duality arguments and deals with 
unbounded coagulation kernels. 

All the work discussed so far deals with diffusing particles (contrast ([2])) and solutions either 
with a zero gradient boundary conditions, which excludes outflow or defined on the whole of so 
that outflow is thereby excluded. For numerical reasons motivated by applications in engineering, 
the present author has been interested in the Smoluchowski equation with advective transport and a 
delocalised coagulation interaction jll, [lcj] • In particular for engineering applications particle gain and 
loss terms are important—industrial equipment is designed to take in material, alter it and then send 
it on either as waste or product. This gives the problem as formulated in 15, 0 and other applied 


works a different structure to those studied in previous mathematical works. For example, individual 
particles experience irreversible processes, but nevertheless the system is expected to reach a steady 
state in the large time limit under a wide range of conditions. Measure valued processes (which can 
be interpreted as particle processes) with an inflow term although no interaction were also studied in 
0 - 

For ([2]) specific problem an initial existence result via the compactness of approximating stochastic 
processes was given in 14]. In that work however convergence of the approximating processes could not 


be proved, only sequential compactness, because the number of distinct limit points was unknown. 
This was not only mathematically frustrating, but also a major obstacle hindering the numerical 
analysis of the associated simulation methods. 

The purpose of the present work is to establish uniqueness of measure valued solutions for (J2]). 
Additionally Lipschitz continuity in the initial conditions is shown and the same Picard iteration 
method that proves uniqueness of solutions provides a purely analytic existence proof. The result can 
thus be characterised as one of “well posedness”. Formally there are some new existence results—the 
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assumption of only one spatial dimension in 14] is relaxed, but with the assumptions used in this work 
the proof in that paper could easily be extended. The existence of a differentiable strong solution is of 
interest, because it opens the way to a study of the way in which the solution approaches a solution 
of the corresponding equation with a local coagulation interaction, see for example 


m 


2 Statement of Main Results 

In order to make a precise statement it is first necessary to go into details regarding the various objects 
appearing in ([2]). The basic spaces are the particle position and type spaces X and y respectively. The 
type space, which carries information about the mass and any other internal details of a particle is 
assumed to be a locally compact, second countable Hausdorff space on which coagulation is represented 
by a commutative + operator. The particle position space X is assumed to be a simply connected, 
relatively compact subset of M. d , which is equipped with Lebesgue measure and a derivative V. Both 
X and y are given their respective Borel cr-algebras and X x y is given the product topology and 
cr-algebra. 

Throughout this work M. d will be given the usual Euclidean norm, which will be written |-|. Linear 
operators L between two normed spaces (A, ||-|| A ) ( B , ||-|| s ) are given the operator norm 

\\ l \\a->b : = SU P W Lx W B- 

xGA: ||a;|| A =l 


2.1 Properties of the Flow and Spatial Domain 
Velocity field 

Particles are assume to be transported in a time dependent velocity field ut defined on X the closure 
of X such that u £ C (M + ,C 2 (T,satisfying 


Mloo : = su PieR+,xeT M x )l < 0°. 


IIV • ulloo := sup teR + xe7P £fc=i Qr k u k , t (x) < oo, 

|||Vu||| := supj ||< oo viewing the matrices Vut as linear operators, 
IIVV • uWn := sup igR+ x |V (V • ut(®))| < oo. 


Boundaries 

It is assumed that the spatial domain X is simply connected and has a regular boundary dX that can 
be decomposed into three parts, each with outward normal n(x): 

• Tin where n(x ) • ut(x) < 0 for all t £ M + , 

• Tside where n(x) ■ ut(x) = 0, 

• T out where n(x) ■ ut(x ) > 0. 

Further T m C X but T s ide, Tout C \ X . 


Flow Field 

Define as the position at time t of a particle moving with the velocity field u starting from x 

at time s. It is assumed that 

• There exists a to > 0 such that, for all f > 0 and x £ X one has &t,t+t 0 ( x ) ^ X, that is, an upper 
bounded on the residence time. 

• For every t > 0 and x £ X there exist unique s(t,x), £(t,x) such that $ s u x ),t (^(^ a: )) = x an< ^ 
either s(t,x ) = 0 or £(t,x) £ Ti n (the possibility of both is not excluded). This defines a start 
position for each point in the flow and £(t,x) = &t,s(x)( x )- 
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• s(t,x) and f(t,x) are differentiable in x and HVs]]^ := sup £a , |Vs(/,x)| < oo. A bound for the 
derivative of £ is given in the appendix. 

• The set E t = {x £ X : £(t, x) £ Ti n } forms a differentiable d— 1 dimensional manifold that divides 
X \ Et into two disjoint simply connected components. 

2.2 Test Function Spaces 

Definition 1. Let B\fy) be the space of bounded measurable functions on y with the supremum ( not 
essential supremum) norm, which will be written |HI;y—oo - 

Definition 2. Let B := B]fX x y) be the space of bounded measurable functions on X x y with the 
supremum (not essential supremum) norm, which will be written ||-|| B . 

Definition 3. Let B d be the space of d-dimensional vector valued functions with components in B. 
This will be given the norm ||/|| B d := sup X)J/ \f(x,y)\, where |-| is the Euclidean norm on W l . 

To handle the derivative in ((2D and associated boundary condition introduce 

Definition 4. 


D := < / 6 B : f differentiable , V/ £ B c 


. ||/(*,-)|| 

x — >r out 


y~oo 


= 0 


The norm is 

ll/llu = ll/llfl + ||V/|| fl<J . 


This is an appropriate class of test functions to use in ©, because the derivative is well behaved. 
For a discussion of the boundary condition see 14], although that work imposes slightly stricter 
regularity conditions, which are here seen to be unnecessary. 


2.3 Solution Spaces 

A particle distribution is at a minimum a measure on the product of the particle position and type 
spaces, that is on X x y. The solution processes must accordingly take values in the following spaces, 
which are built from the space of measures on particle types T 

Definition 5. Let (Af(T’), ll'lly-Tv) =: A4 (T)tv be the normed space of signed bounded measures on 
y with the total variation norm 


W/J-Wy-TV : = SU P 

/AO 


J y f(y)v( d y) 

\\f\\y-oo 


f e B h (y). 


Definition 6. Let M = A4(X x T) be the vector space of bounded signed measures on X x y. 

Under reasonable assumptions one expects to find solutions to ([2]) that are absolutely continuous 
with respect to Lebesgue measure on X; this leads to the following space (compare 01). 

Definition 7. 


M 0 ,oo = {p £ M : p(dx, d y) = c(x, dy)dx, c £ L°° (X, A4(T)tv)} 

with the norm 

ll c IUooc = ®ss sup ||c(x, •)||^;_tv 

X 

where a measure is identified with its density. 

The B and D dual norms on A4 will play a role in this work 
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Definition 8. Let /j, G M. 


TV — 


B~> 


:= sup 
J¥0 


fx x yf( x 'y)v( dx ' d y) 


B 


f£B, 


and 


11/4D* : = SU P 

m 


fxxyf( x ’y)v( dx ’ d y) 

UV D 


f€D. 


As the notation suggests, the B* norm is the total variation norm on At. For calculations the B* 
point of view is emphasised, however the main results are stated in terms of TV. When dealing with 
processes the following abbreviation is useful 


Definition 9. Let T > 0 and c G L°° ([0, T); At) then 


B* 


ess sup ||c(f, -, -)|| B * ■ 
te[o,T) 


2.4 Coagulation 


It is now possible to set out the assumptions on the coagulation dynamics specified by K and h 
in ©• K is assumed to be non-negative and measurable with some bound K 0Q > 0 such that 
su Pt/i,?/2 K (yi>y2) < Koo. 

The delocalisation h: X 2 M must be measurable and non-negative. For fixed x\ G X write h\ xi 
and /i 2 ,si f° r Hre the functions given by hi tXl (-) = h(x i,-) and h 2 , Xl (') = h(-,xi). It will be assumed 
that neither K nor h are identically zero—this would lead to a trivial problem with no coagulation. 

HI: \\h itX \\ B < Ci \/x G X, i = 1,2. 

H2: HI holds and h(x,x 2 ) = YSj=i Xj,i( x )Xj, 2 (% 2 ) with Xj,i positive, and of special bounded 
variation (derivative in L 1 plus atoms) for all i and j with the number of atoms in the weak derivatives 
bounded. Further one has sup XjX2 Ylj=i IIXj,2|ls ft |Vxj,i ($ r , s (a;), x 2 )\ ds < mm (*-r,to) an d 

B St |VXj ,2 (* 1 , 

%hi, x (S) 

here). 


,$ r , s (a;))| ds < C 2 f 0 e lliV “ ll|min(t ~ rA) 

< C 2 Vx G X,i = 1,2. It should be noted 


its symmetric counterpart sup :1 , ix )Cj=i ||Xj,i 
H3: HI holds and the h l;X are in D with 

that H3 implies H2 (Proposition [60] is helpfu’ 

The function h parametrises the numerical methods that lie behind this work [15j]. H2 is describes 
the case where the spatial domain is partitioned into cells and coagulation is only simulated between 
particles that are in the same cell. From a software point of view this is somewhat simpler than dealing 
with functions satisfying H3. In one dimension, which was the case simulated in [15], H2 is a weak 
integrability condition on the derivative of h. 


2.5 Inception 

Particles are added to the system with intensity given by signed measures It G At. 

II: sup t ||/i|| B * < 00 and I G C([0, 00 ), (At, IHId*))- 

12: II holds, the It are non-negative measures and for every / G B 

/ f(x,y)I t (dx,dy) = f(x,y)I mt (t,x,dy)dx + /(£, y) Ibdry (t, C dy)d£ 

Jxxy Jxxy ^r in xv 

with /i nt G C ([0, 00 ), A4o,oo) also Ibdry G C ([0, 00 ), L°° (Tin, A4 (V)tv)) with the respective norms 
uniformly bounded for all time and with some I* > 0 such that 11Ibdry(^j £1 Olly-TV — I*Ut(£) ■ n(£) 
for all t and £ G Tin- 

13: 12 holds, Ibdry has a time and space derivative so that Ibdry € C 1 ([0, 00 ) x ri n , A4 (V)tv) and 
lint has an T-derivative which is VIi nt G C ([0, 00 ), (Ato >0 o) d ) 


5 














2.6 Statements of the Theorems 


These results progress from local existence and uniqueness of a measure valued solution to a global 
result and then existence followed by differentiability of a density for the measures. 

Theorem 10. Assume H2 or H3 holds and that co £ M, then there exists a T = T(c o) such there is 
a unique solution ct to ([2]) in L°° ([0, T), (M, ||• || TV )) with initial condition cq and this solution is in 
C ([0, T), ( M , II-IItv)) n C 1 ((0, T), (M, ||-|| D *))- 

Additionally, there is no time interval on which more than one TV-hounded solution exists for a 
given initial condition. If solutions exist on a common compact time interval for at two or more initial 
conditions, then the solutions are Lipschitz continuous with respect to the initial data in the TV-norm 
on this compact time interval. 

In the physically reasonable setting of non-negative particle numbers, the previous result holds for 
all time: 

Theorem 11. The T = T(co) from the previous theorem is oo if Co and the It are non-negative 
measures. 

Theorem 12. Assume H2 or H3 holds, that cq is in the positive cone of Mq ) 00 , and that 12 is 
satisfied, then fl5J) has a unique solution, which is in L°° ([0, oo), -Mo,oo) and therefore has a density in 
L°° ([0, oo) x X,M.(y) tv) starting from co- 

Theorem 13. Assume that cq £ W 1,QO (X,M.(y) tv) is consistent with the boundary condition given 
below, that 13 is satisfied, and further that either H3 holds and X has a sufficiently regular bound¬ 
ary or d = 1, H2 holds and u is bounded away from 0, then © has a unique solution c with a 
density in W 1,oc ([0, oo) x X,M(y) tv); satisfying the boundary condition —ut(x ) • n(x)c(t, x, dy) = 
Ibdry(t,x,dy) Vf £ M + , a: £ Tin and with initial condition co- 

As a corollary of the preceding two results an earlier result by the author, which demonstrated the 
existence of converging sub-sequences of stochastic approximations to solutions (J2]) can be extended 
to a full convergence result: 

Theorem 14. The stochastic jump processes studied in a which have X = [0, L ) for some L > 0 
and satisfy H2 and 13 converge to the unique solution of © and this weak solution is also a strong 
solution in the Sobolev space IT 1 ’ 00 ([0, oo) x [0, L), A4(T)tv) provided that the initial condition co 
is in IT 1,00 ([0, L), -M(T)tv) with rto(0)co(0, dy) = Zbdry(0, 0, dy). Further one has ut(0)c(t,0,dy) = 
4dr y (i,0,dy) for all t. 

Proof. In 0 it was shown that every sequence of approximating processes has a sub-sequence con¬ 
verging to a solution of © (a compactness result). Theorem 1111 shows that there is only one such 
limit point so one has convergence and Theorem 1131 yields the differentiability. □ 


3 Dual Operator Estimates 


Introduce the more compact notation (/, y) = J Xx y /(x, y)y(dx, dy) for / £ B and y £ M.. It is now 
helpful to seek a generator for the evolution given in ([2]), that is an operator At such that 


_d 
d t 


(f,ta t ) = (Mf),iit) + (f,It) 


(3) 


This is in fact a dual generator, because it acts on the functions not the measures. 

The author emphasises his dependence on Kolokoltsov [j] for the material in this section and the 
first half of the next. The first novelty in this section is the boundary condition associated with the 
finite domain and outflow, which required careful treatment, but is not covered by the existing work. 
Also the consideration of coefficients of bounded variation (H2) is essential to treating the motivating 
example from [1J] and even under this relatively weak assumption differentiability of the solutions in 
one spatial dimension is established. An additional variation from [7] appears in Proposition 1261 where 
some additional problem structure is exploited and enables the fixed point methods to be applied in 
the £>*-norrri; rather than the weaker ZT-norm used in @]. 
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3.1 The Generators 


Because © is quadratic in p the same must be true of the expression (A t [p](f),p t ), which is achieved 
by including the path (p r )re[o,t] as a parameter of A t . It is technically convenient to parametrise by 
the entire path, not just pt, because one eventually deals with propagators where the dependence 
cannot be expressed in terms of p at any finite set of time points. One notes that A t [p] = JJt + H t [p] 
where U is the transport operator and H is the coagulation operator. 

Definition 15. Let p £ L°° ([0,T), (M, ||-|| B *)) and assume HI holds. The coagulation generator 
parametrised by p is H t [p]: B —>■ B defined by 

H t[p]{f)(x,y) = - I h(x,x 2 )f(x,y+ y 2 )K(y,y 2 )p t (dx 2 ,dy 2 ) 

A Jxxy 

- x / f{x,y)[h(x,x 2 ) + h{x 2 ,x)\K(y,y 2 )p t (dx 2 ,dy 2 ) (4) 

1 Jxxy 

fort € [0,T). 

This is not the only possible definition for Ht\p\, other versions also yield the desired expression 
(the coagulation term from ([2])) for (H t [p\(f),p t ). Each definition would lead to characterising the 
solutions as fixed points of a different mapping; the definition given here seems to be the one that 
minimises the technical difficulties in the following analysis. 

Proposition 16. Let 0 < T and p € L°° ([0, T), A4) and assume HI holds, then the operator norm of 
H t [p] as a mapping B —>• B satisfies 

3 

ess sup 11 Hf [p] 11< -ITooCiIIHIIb*- 
t 2 

Proof. Immediate. □ 

Definition 17. Let t £ M and f £ D then the transport generator Ut- D —> B is given by 

U t f(x,y ) = u t (x ) • Vf(x,y). 

One can now define Afip] = Ut + Hfip] as a linear operator D —>• B. 


3.2 The Propagators 

Propagators are generalisations of semi-groups to deal with time dependent generators. For a detailed 
discussion the reader is referred to [7, Chapter 2] or 17, Chapter 5]. The key idea (given in the dual 
setting appropriate to this section) is that a generator At generates a family of linear operators A r,s 
such that A r,s A s,t = A r,t and 


— A^ = A s,t A t , — ,4 s ’* = -AsA^. 
d t ds 


(5) 


The goal of this section is to construct such a family of propagators for the generator Afip] from the 
previous section. 

Definition 18. Let s,t £ R and f € B and define the transport propagators U t,s : B —>• B by 

otherwise 

where T is the flow due to the velocity field u (see { i2.1\ and Appendix^). 

Proposition 19. Let t > s, then the transport propagator U s,t preserves D. The following operator 
norm estimates hold: 

\\ uS,t \\ B ^ B < t(t-s <t 0 ), 
where 1 is an indicator function and 

||P s,t || Z) ^ D < e^- s )IH Vu llll (t-s<t 0 ). 
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Proof. The D-norrn of / is immediate. Use the chain rule and Proposition [60] in the appendix for the 
derivative of /. □ 

The required propagator is now constructed as a perturbation of the transport propagator U by 
the bounded coagulation generator: 

Definition 20. Let T>0,0<r<t<T and p <E L°° ([0,T), M), and define (compare ]% Theorem 
2.9]): 


A r ^[p\ := U r >* + Y, / U r ’ s 'H sl [p)..-U s ™-'’ s ™H Sm [p\U s ™’ t ds 1 ---ds m . 

It is now necessary to establish estimates for the operator norm of A on B and D. For this it is 
shown that the infinite sum just given is absolutely convergent in both operator norms. During this 
analysis it is convenient to use some additional notation: 

Definition 21. Under the assumptions of Definitional^ let f € B and t > 0; define both f® t := U r,t f 
and 

f™ t := 

J r 

This allows one to write 

OO 

A r, w = E /”«• m 

m =0 

Proposition 22. Under the assumptions of Definitions \WJf ; {27\ 

1/1 \ m 

which is zero for t — r >to and 

1/1 \ m 

Proof. Proceed by induction. □ 

The ^-operator norm estimate now follows: 

Proposition 23. Let T>0, 0 < r < t < T and p, € L°° ([0, T), (M, ||-|| s *)) and assume HI holds, 
then A r,t [p] is a locally bounded propagator on B satisfying 

\\ Ar,t H\ b^b ^ - r <t 0 ) 

and A r ’ t [p]f is \\\\ B -continuous in t for every f G B and t > r (this is known as ‘strong continuity'). 
Further, for any f £ D and almost all t 

= A^\p]A t [p\f, 

where one recalls A t [p\ = ut • V + H t [p\. 

Proof. The first part of the result follows from Proposition [22] and ([6]). 

The (left) generator At[p] can be found differentiating the series in Proposition [20] term by term 
and observing that the resulting series is again absolutely convergent. □ 

Differentiating with respect to r in Proposition [23] is not possible, because A r,t does not necessarily 
preserve D. This is addressed in the next few propositions by making stronger smoothness assumptions 
on h, the spatial delocalisation of the coagulation interaction introduced in 112.41 and used in the 
definition of H (Definition fl5l) . 
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Proposition 24. Under the assumptions of Definitions \2(MA21\ and additionally assuming either H3 
holds or H2 holds and // is bounded for all (not just Lebesgue almost all t), for example because it is 
continuous 


v/ r m t(av)|| y _ oo < ( -ifocAllMll 


I B* 


(t - r) r 


ml 


|V(/(‘M*)>-))lly-c 


'3 \ m (t - r )™- 1 Cotne^ Vu IIK t_r ) 

+ m [ Ul/llll** ) k J 20 ^ - ||/ -)||y_ c 


Cl 


and 


v/ r m *|| Bd < l-KMMW 


B* 


{t - r) r 


ml 


! e ll|Vw|||(*-r) \\Vf\\ Bd t(t-r>t 0 ) 


3 \ m (+ - r \ m ~ l Cotne^ u ^ t - r ' ) 

+ m ( -KooCMs.j 1 J 20 ^ - \\f\\ B 1 (t-r> to). 


Proof. The first inequality is established by induction making use of Proposition [22] for the terms in 
/. The second inequality introduces Proposition [60] to get an estimate for V4> r) t(x) □ 

Proposition 25. Let T > 0, 0 < r < t < T, p € L°° ([0,T), (AA, ||-|| B *)) and either H3 hold or H2 
hold but with (jl bounded for all (not just Lebesgue almost all) times, then A r,t [fi\ is a propagator on D 
and there is a C 3 € M such that 


\\ Ar,t H\ D ^ D ^ e( |l|v “ lll+ i A '°° Cl|l ^l | ls*)( i - r )c , 3 l(t -r< to). 


Further, for any f € D and almost all t 

^A r ’%]f = -AMA^\p\f. 
Proof. From Proposition 1241 one sees that, for / £ D 


|| V (A r ’*[fi]f) || Bd < e (l^ C llllMlll s * + ll|V«|||)d-r) 11 ^ _ r < to ) x 

(l|V/|| B d + (l + ^KMM\ B .(l - rfj ||/|| B ) . (7) 

For the / part of the D-norm use Proposition l23l the first statement of that proposition also established 
the boundary condition for D. Differentiation in r and t is performed term by term in the infinite 
sum from Definition 1201 □ 

These results concerning the dual propagators are concluded by showing Lipschitz continuity in 
the measure valued path parameter. 

Proposition 26. Let T > 0, suppose n,v £ L°° ([0, T), (M, ||-|| B *)) ; 0 < s < t < T and HI holds, 
then 


lA^p) - A s ^\ 


I B->B 


< ^/iooCie 3XooCimaxl 


B* ’ll 


t)P s )i(4 — s < to) ess sup || n r — ty|| B * 
re[s,t] 


dr. 


and 


lA^lp] - A^u] 


I B^rB 

< ^X 00 Ci(t-s)ei AooCimaxi 


< t 0 ) y* 


— Z/, 


r\\B' 


dr. 
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Proof. Write M = max(|||/r||| B *, 


B +) and show by induction that 


| U r ’ Sl H Sl fi] ■ ■ ■ U Sm ~ 1,Sm H Sm [p\U Sm,t - U r,Sl H Sl [v] • • • \\ B ^ B 

m —1 m 


< r^ooC'i -K^CxM 


— v< 


j 11 B* 


3 = 1 


( 8 ) 

□ 


This result exploits a small amount of additional problem structure to adapt the method set out 
in the proof of Theorem 2.12 in [?]. The key is that the parameterisation only affects the coagulation 
( H ) part of the propagator, which has a bounded generator, while the transport ( U) part of the 
propagator, which has an unbounded generator is independent of the parameterisation by p and u. 


4 Operators on the Space of Measures 


Under the duality pairing of B and A4 given by (f,p) = f Xx y fp(dx, dy) as used above, (dual) 
operators B —>• B define (pre-dual) operators AA —> M with the same operator norms. 

Definition 27. Let 0 < s < t < T and /a € L°° ([0,T), (A4, ||-|| B *)). For the pre-duals of U s,t and 
A s,t [p\ write U t,s and A t,s [p\ respectively and note the reversal of the time indices. For the pre-dual 
of H t [p] write H t [p]. 

It is emphasised that A s,t [p\ acts on functions while A t,s [p] acts on measures, but both are param- 
eterised by a measure-valued path p. 

The existence of the dual operators and their norm estimates is immediate, see for example [7, 
Thrm 2.10]. The duality relations yield: 

Proposition 28. Let 0 < s < t < T, p,v £ L°° ([0, T), M) and assume HI holds, then 




M->M 


= \\A‘M\ 


B->B — 


e |ftrooCi|||//||| B *(i-s) 1 (t - S <to). 


along with 


A^lul - A^ s [u] = || A s ^[p} - A s ^[u] || 


M^M 




< ■^K 00 Cie 3KcoCl max (IH^IIlB*.Ill^llls* - s <t 0 ) 


re[s,t] 


\\Pr ~ Vr\\ B * dr ‘ 


and 


A^ s [p] - A^ s [u\ 


M^M 
3 


= ||a s >V] -a^mw 


B-s-B 


< jrKMit - s)e3 Jf ‘» ClI1M x(IWII B *,|ll»'lllfl*)(t-«) 1 (t _ s < t 0 ) ess sup || p r - v, 


r6[s,t] 


r\\B* 


Proof. Duality and Proposition 1231 


□ 


Proposition 29. Let 0 <s<t<T,p£ L°° ([0, T), (A4, ||-|| B *)), c € M., f G D and assume HI 
holds, then for almost all t 

±(f,A t ’ s [p]c) = (A t [p}f,A t ’ s [p]c)^ 

Proof. Duality and Proposition [23] □ 
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4.1 The Fixed Point Mapping 


This section presents a Picard iteration method for @ highlighting the roles of the B* and D* norms 
on the space of measures. The mapping that will be shown to have a fixed point is: 


Definition 30. Suppose Co £ A4, 0 < t < T, let /i £ L°° ([0,T), (Ad, ||-|| s *)) and suppose HI holds. 
Define T Co : L°° ([0,T), (Ad, ||-|| B *)) ->■ L°° ([0,T), (Ad, ||-|| B *)) by 




co 


(vm = A t ’°[ f i]c 0 + f A^Wsds. 
Jo 


Proposition 31. Under the assumptions of Definition \30\ one has 'P CQ (p) £ Cb ([0, T], (Ad, 
with 

" ls Ms* 


< e§ K » Clllw "»*“”«•'»> (llcclls. !((-*<(„) + 


B* 


The time derivative exists for almost all t £ (0, T) with 

MW 


D i 


< WhWo* + ( llulloo + 2^00^1 ||/itlls* ) ll^co (u)(t) II B* 


and if n £ C([0,T), (M , \\-\\ B *)) then T co (/r) £ C 1 ((0,T), (Ad, ||-|| D *)) ■ 


(9) 
I B*)) 


( 10 ) 


Proof. B * boundedness is a consequence of Proposition [28] and continuity follows from the continuity 
in t of A t ’°\/j]. 

For the time derivative differentiate the formula in Definition [30] and use Proposition [29] □ 


Proposition 32. Suppose Co £ Ad, T £ (0, oo), H2 or H3 holds and c: [0,T) —> (Ad,||-|| B *) is a 
bounded solution to (ED with initial condition cq, then c is a fixed point of 


Proof. Suppose c to be a solution of ([2]) and let t £ [0, T), then using duality and Proposition 1251 one 
finds 


d_ 

dr 


f,A t ’ r [c] Cr ) = (f,A t ’ r [c]I, 


( 11 ) 


Integrating over r £ [0, t] completes the result. This (standard) argument can be found, for example, 


m 


17, §5.1], 


□ 


Proposition [22] is the only place where one requires H2 or H3 in the existence and uniqueness 
analysis. This is in order to invoke Proposition 1251 and more fundamentally so that A s,t [/j] preserves 
D ; otherwise one cannot give meaning to 4-A t,r \ij]. Without this result it still follows that the 
mapping T has unique fixed point with all the advertised properties (in particular solving ([2|l ). but 
one cannot rule out the possibility that there are additional (possibly less regular) solutions to ([2]). 
These conclusions are stated more formally in Proposition [35] for which two preparatory results are 
needed. 


Proposition 33. Let cq £ Ad, M £ M + be large enough to satisfy 


M > ||co|| B . 


2sup s \\I S \\ B * 
3/Pqo Ci Af ’ 


define Em = {/r £ Ad: ||/i|| B * < M} and assume HI holds. Then there exists a tm > 0 such that T Co 
preserves L°° ([0,r M ), (E M , INIs*))- 


Proof. Let t'm > 1 be given by 




IIcoIIb* + 


2s "Ps II^IIbA 

3/Poo C, A/ ) 


= M 


( 12 ) 


and suppose /r £ L°° ([0, T), (Em, IMIb*)) f° r some T > 0. Use Definition [301 along with the operator 
norm estimate from Proposition [28] to see that, for t <T 


|| B * 


< c %KaoC\M min(i.to) M 

~ ' r M 


(13) 


and so ||4' Co ( / u)(t)|| s * < A I if min(t, to) < 3 ■ Hence it is sufficient to take tm = an d ^ 

to, the maximum residence time for a particle, satisfies to < 3 “k^CiM then one may take tm = oo. □ 
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Proposition 34. Let cq S M. and Em,tm be as in Proposition E3 and assume HI holds, then there 
is a t' m < tm such that is a contraction on L°° ([0 (Em, ||-|| b *)). 

Proof. Suppose /a and v are in L°° ([0, tm), (Em, ||-|| b *)) / £ B and t £ [0, tm), then by Proposition [26] 


||4/ C0 ( M )(t)-vh C0 (^)(t)|| 


IB* 


< 


llcolls. 1 (t < („) 

ID — 


+ [ A t,s [(i\ — A t,s [v\ ||/J B *l(i-s<i 0 )ds 
Jo S-)-B 

< -K 00 Cite% K ° oClMt ||co|| B * t(t < t 0 ) ess sup \\n r - v r \\ B + 

^ r£[0,t M ) 

+ ^KooCi min(f 2 , f^'j e i K oo c i Mimn (. t , t o) gU p ||/ r |[ s * ess sup || /j r — u, 


re[o ,t M ) 


r\\B* ■ 


Hence for any 0 < r < 1 one can find a t' m < tm such that 

sup ||^co(M 1 )W-^co(^ 2 )WL* < resssup \\/4 - ht\\ B 
te [°’ T M) *e[o P M ) 


(14) 


(15) 

□ 


Proposition 35. Let Co £ A4 and Em be as in Proposition \3SX t' m as in Proposition [3J] and as¬ 
sume HI holds, then (J2[) with initial condition co has a solution on [0, t' m ) and this solution is in 
Co ([ 0 ,t' m ), (Em, ||-|| b *)). If H2 or H3 hold this solution is unique. 

Proof. By Propositionl34l there is precisely one fixed point of ’Peg, which by Proposition [31] is a solution 
of ([2jl with initial condition co- Proposition [33] shows that this solution is in C\, ([0,t^-) ? (Em, IHIb* ))■ 
By Proposition [32] every solution of ([2]) with initial condition co is a Hxed point of 4% and thus is 
unique. □ 

Proposition 36. Let T > 0, assume HI and suppose l I r /i0 and 4 , l/Q have fixed points /i and u respec¬ 
tively. Write M = max (|||/i|[| B +, |||v[|| B *), then there exists C±(M) > 0 such that for t <T 

||ht - Vt\\B* < Who - Ms* 

and thus at most one finite solution is possible for any given initial condition. 

Proof. Since any solution must be a fixed point of 1 P for the appropriate initial condition 


II ht ~ u t\\ B * = ||^ 0 (m)W - ^v 0 (v)(t)\\ B * 

< II^mo (h)(t) ~ ^oHWII B * + - ^oHWII B * • ( 16 ) 

Now by Proposition [28] estimate the second term as follows 


II *«»(*) - ^M(0IIb* 


A tfi M (ho - Vo) 


< e i K °° c i Mt t (t < to) Who 


vo\\B*- (17) 


For the first term using Proposition [28] one finds 


(h)(t) ~ 


IB* 


< 


n „ ll^olls* 1 (t < to) 


ft,Or 


B—»B 


+ [ A*’ 3 [fj] - A*’ 3 [v] ||/ S || B * 1 (t - s < t 0 )ds 

Jo B ^ B 


< — Kqo C\ e 3K °° ClMt 


B’ 


+ B<ooCW 


1 (t < to) [ II hr 
Jr£[0,T) 

SKooCiMmin(t,to) T ni n ( <] gup || ^ 


— v r \\ B t dr 


/ 

Jr€ 


r\\ B * I ||hr 11 B * - (18) 

re[0,T) 
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so using Gronwall with 


C 4 (M) = ^K ao C ie v K °° ClMt0 (||cS|| B * +i 0 sup||/r|| B *) (!9) 

one has 

II fH - Vt\\ B * < WhO - Mb* e I^ClMmin(^o) e C 4 (M)t_ (20 ) 

□ 

Proof of Theorem 1 HA The existence of a solution on a small time interval is the conclusion of Propo¬ 
sition [35l this procedure may be iterated, but the time steps may decay so that a solution cannot 
necessarily be constructed for all time. 

Proposition 1321 establishes a representation for any solutions, should they exist. Using this rep¬ 
resentation boundedness and continuity in the £?*-norm along with differentiability in the D*-norm 
were established in Proposition 1311 

For compact subsets of the time interval on which a solution exists (which may be longer than the 
time interval for which this theorem proves existence), B* Lipschitz continuity in the initial conditions 
and uniqueness are consequences of Proposition 1361 □ 


4.2 Positive Measures 

Write B + for the cone of non-negative functions in B and M + ,A4q 00 for the cone of non-negative 
measures in At, respectively A4o,oo- These cones are of course not Banach spaces, but one would 
expect the physical solutions of any reaction-transport problem to remain in A4 + , if they start there. 
This is indeed the case and turns out to allow the local existence result for the coagulation-transport 
problem studied here to be extended to a global one, which along with the results already established 
makes the problem well posed. 


Proposition 37. Let T > 0 and [v E L°° ([0, T), (A4 + , ||-|| B *)), then for 0 <s<t<T A s,t [/j] is a 
positivity preserving on B, the same is true of A t,s [fj] on M. and both operators are contractions on 
the respective positive cones, that is 


K'WIIb < 


Bi 


/ e B- 


X' a W 


B* 


— IMI B* ) 


v E M + . 


Proof. A proof for the dual propagators on B suffices. For this note that U s,t is positivity preserving 
with P-operator norm 1. One further checks that Ht[p] generates a positivity preserving propagator 
with operator norm at most 1 on B + , which will be denoted H s,t [/j,\. One can now approximate A s,t [fi] 
by 

t . = s + i Lll j i = l,...m-l, men (21) 


which is a splitting, to see positivity is preserved and the operator norm is bounded above by 1. □ 


The key estimate from Proposition ED can now be improved (recall to is the maximum particle 
residence time from H2.ll) : 

Proposition 38. Assume HI holds, cq E A4 + and /x G L°° ([0, T), A4 + ) for T G [0, oo) then 


II ^c 0 (h)(f) II b* < IMIs*l(t < to) + min(f,fo)sup||/ s || B * . 

S 

Proof. This follows from Definition l30l and the norm estimates in Proposition [371 EH 

Proof of Theorem , [771 One can take M = ||co|| B * + tosup s ||/ S || B * and tjw = oo in Proposition [33j 
Proposition [M] then extends to show that T Cq is a contraction on L°° ([0, oo), Em H A4 + ). □ 
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4.3 Measures with Lebesgue Densities 

One would of course like to prove that every measure valued to solution to © is in fact also a strong 
solution to an appropriate extension of ©• The main difficulty that has to be addressed in this 
section is the inflow of pre-existing particles through Ti n which leads to It having a singular (with 
respect to Lebesgue measure on X) part concentrated on Ti n . In this section it is shown that under 
a mild time-regularity condition (12) the advective transport smooths out the inception concentrated 
on Tin sufficiently for solutions to ([2]) to remain in _A4o i0 o- Shocks are of course preserved by advective 
transport, but what happens here is more like spraying paint onto a moving surface, as long as the 
surface keeps moving a thin layer of paint is deposited everywhere and no ridge (shock) is created. 

Proposition 39. Assume HI holds, 0 < s < t < T, //£ L°° ([0, T), (_A/[, ||-|| B *)), then 

oo 

= U l ’ s + ^ ■ ■ ■ H Sl \jj]U r ’ Sl dsi • • • ds m . 

Proof. For each ?n the term in the sum here is dual to the term with the same m. in Definition [20j □ 

Proposition 40. Assume HI holds, 0<s<t<T,y€ L°° ([0, T), (M, IHIb*)) and c € M.q jOQ , then 
for any cf G £>b(T) and bounded measurable f: X —» R 

[ f(x) [ (f(y)U t ’ s c(x,dy)dx = [ f(x)e~ v -“r(*t,r(aO)dr f <^(y)c (4> M (x), dy) dx 
Jx Jy Jx Jy 

and 


i.: 


4>(y) H t[lAc(x,dy) = 


ly z Jxxy 


(t>{y + V 2 )h(x, x 2 )K(y , y 2 )nt(dx 2 ,dy 2 )c(x, dy) 


- x/ <t>(y)[h(x,x 2 )+ h(x 2 ,x)\K(y,y 2 )y t (dx 2 ,dy 2 )c{x,dy). 

Jy z Jxxy 

Proof. For the first statement, which concerns the transport propagator U, one makes the change of 
variable x -f->- 4> Si t(x). Liouville’s formula then gives the determinant of the Jacobian as det = 

exp V • ut (4 , t !r (x)) dr. Alternatively one can approximate c by A-differentiable functions (since the 
claim is only of an L 1 nature) and check the formula directly using Proposition 1591 

For Hf use Definition [T2I the important point is that the new measure also has a density with 
respect to Lebesgue measure on X. □ 


Proposition 41. Assume HI holds, 0 <s<t<T, 
bounded propagator on Mo,oo with 


A^[p\ 


-^0,00—^A^ 0 ,oo 

Proof. From Proposition 0D] one sees that 


0, T), (J\4, |HIb*)) then A t,s [y ] is a 

< e (llV-u|| 00 +fx 00 C 1 ||| / ,||| s *)(t- s ) 1 ( t _ < tQ y 


u 


t,s 


Mo , oo —*Mo,o 


< e 


IV'UlUt-8) 


and 


HtilA 


< ^kmmWb* 

-^1o,oo—^A^0,oo " 


The proof now follows that of Proposition [23j 


□ 


The U and therefore also the A are not (norm-)continuous on _Mo,oo- This can easily be seen by 
considering a small translation of a step function regarded as the density of a measure in A4o,oo- The 
eventual time continuity of the solutions will depend on having some T-regularity for the densities of 
the measures in A4o,oo- 

The next proposition provides a better norm estimate when the propagator is restricted to positive 
measures. This is then used in Propositions FiTilD l ll to show that inception concentrated on the inflow 
boundary does not take the solution out of Mq,oo- 
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Proposition 42. Assume HI holds, 0<s<t<T,/j.£ L°° ([0, T), (_A4 + , ||-|| s *)) and let c £ A4o,< 
be a positive measure, then A t,s [/j]c is also a positive measure and 


A*’>]c < e llv-IU(t-) 1(i - s <to) ||c| 


Mo. 


Mo :00 


Proof. Since Mo,oo C M preservation of positivity is a consequence of Proposition [37] and Proposi- 
tion0O]states the A4o,oo is preserved. To proceed note that the propagator generated by H t [p\ preserves 
Mq 00 not just A4 q,oo and coagulation reduces c(x, y) = \\c(x, -)||y_ T v for all x £ X. Secondly U t,s 


is positivity preserving and 


U ] 


t,s 


Mo,a 


< using the representation from Proposition 131)1 


so 


the result now follows by the same splitting approximation as in the proof of Proposition [371 


□ 


Proposition 43. Let T > 0, and v. [0, T) M be such that (note the reduction in the domain of 
integration accompanied by a change in the position of the time argument) 

[ f(x,y)vt(dx,dy)= [ /(f,2/Mi,£,dy)d£ V/ £ B Vfe[0,T). 

Jxxv ./r„,xv 


Suppose further that there is a v* £ (0, oo) such that sup ig [ 0jT ),^gr in ||^(^£j Olly-TV / u t{0 ' n (0 < u *> 
then Jq U t,s n s ds £ L°° ([0, T), A4o,oo) and for all (not just almost all) t <T 


f 

Jo 




< V e ll V " Lt lloo m in(i,t°) _ 


Mo,a 


Proof. Let ( £ P m and take an orthonormal basis for M. d at £ given by e\ = n(£) the outward 
normal and ,ed £ Tin- With respect to this basis let the rows of the matrix V^ r ,t(x) | x= £ be 

diA> r j{x) |a; = g. Thus rows 2 ,d of this matrix are the same as rows 2 ,d of and using 

Proposition [33] the hrst row is 

d 

= V‘h r - i t(x) [x=f ^r(C) = 

d d 

^ ^ ^ \x =£ (^2 * (0) = ^ ^ * ^r(O) » (22) 


2—1 


2=1 


which is ± (di<§> ri t(x) U=<j) ( n (£) * u r(€)) pl us a li near combination of the remaining rows. One thus 
has for £ G Tin 

d®r,t(Q\ = _ Ur (Q . n (£)det (V$ r ,t(x)) 


det 


d(r,£) 


(23) 


x=£ 


Now let f £ B with / ($ r ,t(€),y) = 0 for ^ X as in the definition of U r,t so 


f, f U t,r v r dr\ = f [ f f (<f> r ,t(.0,y)v(r,t,dy)dZ<h 
Js / Js J r in Jy 


Jx: s=$r,t(0 
r€(s,t),i& r in 


det 


d*rM) 


-1 


f (x,y)v(r,£,dy)dx 


Jx: X=<b r ,t(£) 
re(s,t),$er in 


d{r,0 J Jy 

det (\7d>r t t{f)y l 11/ (x, ■)lly_oo K(£) • ^(/)T J Oil dx 


< e" 


IIV-uH^ min(t-s,t 0 ) 


Jx : X=$r,t(£) 

r£(s,t),£e r in 


11/ (x, Olly-oo dx, (24) 


where u(x,dy) is defined to be v{r, £,dy) for the unique r, £ such that & r ,t(£) = x. Proposition l60l in 
the Appendix provides the estimate for the determinant. □ 


15 




















/' 


< ZA e l |V '“ll~ min(t ’ to) . 


Proposition 44. Let T > 0 and ji G L°° ([0,T), (_A/J + , IMIs*)), t/ien under the conditions of Proposi- 
tion\j3\ 

r t _ 

A t ' s [g]n s ds 

Mo,oo~ 

Proof. Use the series expansion from Proposition [39] and the -Mo,oo-operator norm estimates from 
Proposition [42] □ 

Proof of Theorem \1A Theorem |TT] provides the existence of a solution c. Proposition [32] shows that 
this solution satisfies 


c t 


= A t >°[c]c 0 + / A^lcjlsds t> 0 . 


(25) 


Propositions fnhfTTl show that this is in _A/Jo,oo for all times. The boundedness follows from the 
estimates in the same two propositions. □ 

In order to obtain a strong solution to ([2]) it is not sufficient that the measure valued solutions 
have a density with respect to Lebesgue measure on X, this density should itself have a derivative. 


4.4 Differentiability 

One could proceed as in Proposition 1251 to see that Gl t,s [c] preserves measures with T’-differentiable 
densities except for possible jumps where s,t, x are such that 4>t iS (x) G Ti n . This leaves two questions 
open—how to handle these jumps and secondly the treatment of the integral term from f[25l) and 
in particular the /bdry P ar t of I in that integral. The right approach to these tasks seems to be to 
introduce the space of measures with df-bounded variation densities: 

Definition 45. 

Mbv = {cG >fo,oo: (3 C = C{c)) (v/ G D d ) ((V -f,c)<C \\f\\ Bd )} . 

This is equivalent to the existence of a measure Vc G A4 d (not necessarily in Atooo) such that 
(V/, c) = — (/, Vc). 

Until now the notation (/, fi) has been used for f Xx y f(x, y)fj,(dx, dy) for / G B and [i G A4. To 
consider derivatives it is necessary to move to vector valued functions and measures; to facilitate this 
the notation is extended so that for g G B d and v G M d 

d /» 

(g,v)'■= 9 i{x,y)vi(dx,dy). (26) 

i=i Jx *y 

Some more definitions are now needed for the proof that J)* A t,s I s ds and by extension the entire 
solution is -A4 bv- First recall s(t, x) from iJ2.ll the time at which a particle travelling with the flow 
must have entered the domain in order to reach x at time t. Since t is fixed in the relevant places s(x ) 
will be written for brevity in numerous sub- and superscripts, the t should be understood. 

Definition 46. Let f G B, t > 0 and g G C ([0, t], M). Define f^ t = U r,t f, and /^" +1 = 
jt u r ’ s H s [g]f™ t ds as in Definitions^ The define f by 

f™t( x >y) = f™t 1 ($t,r(s) e X) 

the operators S r,t : B —^ B by 

OO 

S r ’ t f(x,y)='£f™ t (x,y) 

m =0 

and finally the operator S l : B —> B by 

OO 

tff&y) = S^f(x, y ) = Y, fstx)A x ’ V)- 

m =0 

This operator can also be regarded as acting on B d by applying it componentwise. 
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Proposition 47. Let f £ B or f £ B d , t > 0 and y £ C ([0,i], M.), then 

l|S‘/(av)|| J ,_ 00 < ||/(x,.)|| y _ 0o 

for all x £ X and the operator even preserves D. 

Proof. This is an exercise in estimating the terms of the summation as in Propositions [2211 1241 □ 

Proposition 48. Let F £ D d , t > 0 and y £ C ([0, t], M), then the operator (V • <S f ) : D d —>• B 
defined by 

(V • S l ) F = V • (5*F) - S* (V • F) 

satisfies 

II (v-5*) < 11^(1, Olly-oo X 

(3^(72(0111^11^. + HVsIU) (l + \k x C, W| b . if - i(j))j 

for all x £ X. Of course S l depends on y, but since this will always be the unique solution to ([2]) this 
detail is ignored in the notation. 

Proof. Define F™ by replacing / with F throughout Definition 0U] and let / = V ■ F and let ff n t be as 
in Definition S6l By induction one establishes 


V-f%(x,-)-j™(x, 


y-oo 


< \\F(X, Oily—oo X 

^K oo C 2 t 0 \\\iJ-\\\ B *m ( ^ooCatolllMlII 


I B* 


m ~ 1 (t - six))™- 1 
(m — 1)! 


• (27) 


One then establishes a similar formula with r replaced by s(x) and the result follows. 


□ 


Proposition 49. Assume H2 or H3 holds, that cq is in the positive cone of Msw, and that 13 is 
satisfied, then the unique solution c to m given by Theorem, UTl satisfies ct € AJbv Vi € M + . 

Proof. Following (1251) and H2.5I c t can be decomposed as 

c t = A t ' 0 [c]c 0 + f A t ’ s [c]I int (s)ds+ f l 4 ’ s [c]/ bdry (s)ds. (28) 

Jo Jo 

One checks in the same way as for the dual propagator in Proposition \2M that 2l f ’ s [c] preserves A4 bv, 
possibly introducing a new jump on the manifold {x £ X : 4>£ )S (x) € Tin}. This deals with the first 
two terms in the above representation; the third term is somewhat more challenging. 

Let F £ D d and t £ M + , then it is sufficient to show that 

f (v-F,A t ’ s [c\I hdiy (s))ds= ft [ {A s ’ t [c}f)^y)L hdry (s,^dy)dCds (29) 

J o N ' Jo J r in Jy 

is bounded by a constant times ||F|| S d. One can introduce a change of variables (s,£) £4 x where (t 
is fixed) $ s ,t(€) = x, so f = 4>t iS (x) is the point where fluid reaching x at time t entered the domain 
and the time of entry was s. The determinant of the Jacobian for this transformation is the inverse of 

det = - e - £ v-u,(**,,(*))( 30 ) 

o{s,0 

by Proposition [60] and the additional factor of —u s (f) ■ n(f) comes from replacing the X direction 
perpendicular to Tin with s (£ lives in the d — 1 dimensional manifold Tin) so 


J*(v-F,A t ’ s [c]I hdTy (s))ds = 

(5‘V-F) (x,y)e^ v ' Ur ^^ x))Ar (u s{x) (^{x)) ■ n(f(®)))J b diy(a(®), £(*), dy)dx 

= (<S*V-F> t ) (31) 



x Jy 
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where 


vt(x,dy) = e^ v ' Ur ^ t ’^ x))dr (u a(a! )(f(®)) ■ n(f(x))) 1 4dry(s(x), f(x), dy). 


(32) 


Writing 

S f V ■ F = V • (5*F) - (V • S 1 ) F (33) 

and checking that u t E A4bv concludes the proof. □ 

To simplify the remainder of this section it will be assumed that X = [0, L ) x Ti n for some L > 0, 
that is, that X is a something rather like a cylinder. The results are expected to generalise, but 
this assumption avoids introducing technical conditions on X. In particular x E X can be written as 
(xi,x 2 ,--- ,x d ) for xi E [0, L) and (x 2 ,--- ,x d ) E Fin- 

Proposition 50. Assume HI and 13 hold, c E L°° ([0, oo), AAo,oc) solves m with additionally c E 
W 1,0 ° ([0, oo) x X,AA(y) tv). Let n{x ) be the outward normal on T m , then 

ui,t(x)c(t, x, dy) = -u t (x) ■ n(x)c(t , x, dy) = Ibdry(C x, dy) Vf E M + , x E T in . 

Proof. Consider fl2J) with f(x i, x 2 ,..., x<j) = {xi < e} as e —> 0. □ 

The existence of a one or more inverses to the divergence operator is necessary to avoid making 
statements about an empty set of functions in the remainder of this section. 

Proposition 51. Let f E D, then there exists g E D d such that V • g = f and g ■ n = 0 on r s id e , 
where n is the outward normal. 

Proof. Take g(xi , x 2 , ■ ■ ■, x d , y) = (gi , £J 2 , - - ■, go) where 

gi(xi,x 2 ,...,x d ,y) = - f(Z,x 2 ,...,x d ,y)d£ (34) 

J X\ 

and gi = 0 for i > 1. This construction has a natural generalisation in terms of path integrals. It is 
not important exactly which end point on T out is chosen because / = 0 all along this boundary. □ 

This representation is not in general unique. Consider for example the case where / = 0 on T^je 
and take integrals along lines perpendicular to the direction used in the above proof. 

If c solves d2J) , / Efl and g: X —> C\,{y) d is differentiable with V ■ g = f and g ■ n = 0 on T si( j e 
for normal vectors n, then applying the divergence theorem to ((2|) (at this stage in a purely formal 
calculation) suggests 


— f g-Vcdxdy= f uj (V g)V c dxdy — f g T (V ut)V c dxdy 

dt Jxxy Jxxy Jxxy 

- / g • V(V • u t )cdxdy + / g-X7Ii nt dxdy 
JXxy Jxxy 


+ 


L 


r in xy 


g ■ n^-cdxdy - 
dt 


uj (Vg)ncdxdy — 


g ■ nI [nf dxdy 


L 


xy 


J FinXy 

v • ff/bdry dxdy + 


T in xy 


g T (Vut)nc dxdy 


H in xy 


+ 


ixxyJr in xy 


K{y , y 2 )y(x, y + y 2 ) ■ n(x) (h(x, x 2 )q(x, dy)) c t (x 2 , dy 2 )dx 2 dx 


- x / / K(y,y 2 )g(x,y + y 2 ) • V (h(x, x 2 )c t (x, dy)) c t (x 2 , dy 2 )dx 2 dx 

z JxxyJ Xxy 

l [ [ K(y,y 2 )g(x,y) ■ n(x) {{h(x,x 2 ) + h(x 2 ,x)) ct(x,dy)) ct(x 2 ,dy 2 )dx 2 dx 

1 JxxyJr in xy 

/ /v(y,y 2 )y(x,y) • V ((/?.(x,x 2 ) + h(x 2 ,x)) c t (x, dy)) c t (x 2 , dy 2 )dx 2 dx. (35) 


1 

+ 2 


IxxyJxxy 


Definition 52. Define a norm on M d , which by a slight abuse of notation will also be referred to as 
the B*-norm by setting ||/i|| B * = sup /eBd . ||/|| fl =i \(f,v)\, for g E M d . 
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Conditions are now provided to make ([351) rigorous, first by restricting the test functions to the 
interior of the domain, so that the boundary terms can be ignored and then proceeding to more general 
test functions: 

Proposition 53. Assume H2 and 13 hold and that c £ L°° ([0, oo), Ado,oo) solves (j2|). Suppose further 
that ct £ -Mbv for each t, so that there exist vector measures u t of finite total variation such that 
(V • /, cf) = — (/, vt) for all f £ D d . This means that v £ L°° ([0, oo), (X4 d , ||-|| B * )) and in particular 
for all g £ (X°, £>b(T ) d ), the space of once continuously differentiable functions with compact 

support strictly contained in the interior X° of X that also satisfy V • g £ D (for example g £ 
cl{x°,B h (y) d )) one has 



(( ^g) T ut,vt ) - ((Vtt t ) T ssz^ - (g ■ V(V • u) ,ct) + (.s, V/ int ,t) + (H t {c]g, 

+ 77/ I g(x,y+ y 2 )K(y,y 2 )ct(x,dy) ■ (V x h(x^)(dx))c t (^dy 2 )df, 
z JxxyJxxy 

~ \ [ [ g(x,y)K(y,y 2 )c t (x,dy) ■ (V x [h(x,£)+ h(£,x)\(dx))ct(£,dy 2 )d£. 

1 JxxyJxxy 


Here H t [c] is a bounded linear operator mapping B —>• B acting componentwise on g with 

an d Ht[c\g(x,y) = H t [c\g(x,y) (recall H t is specified in Dehnition \T5f) for all y and all x 
except possibly x at which h and c both have discontinuities, which is a set of (R d -Lebesgue) measure 
0. Similarly ||c||.m 0oo = ll c II.Mooo w possible differences between c andc on the same set of measure 
0. Because h is only assumed to be of bounded variation, it only has a weak derivative; in the case of 
the weak derivative with respect to the first argument this is written V X h(x,f)(dx). 



Proof. The boundary integrals on T vanish because g is zero here. Note that the product of two 
functions of bounded variation (c and the h in the definition of H) is itself of bounded variation, 
but the Leibniz product rule for differentiation has to be adapted slightly at points where both are 
discontinuous (yielding H and c). The details follow from [2|, Theorem 3.96 & Example 3.97]. In one 
dimension this amounts to adjustments to give left or right continuity at the jump points. □ 


The terms in the preceding expression can be grouped as follows (c is in this context known): 

• Transport ( uJVg,ut ) = (Utg,v t ) (note \7g is a matrix), 

• linear reactions — ( gVut , vf) + (^H t [c\g, utj, 

• source terms, which are collected as a vector measure Jt[c\ so that, for g £ Ck (X°, B\fiy) d ) 


(< 9,Jt[c \) = - {g ■ V(V • u) ,c t ) + (g, V/int.t) 

+ 77 / [ g( x ,y + y 2 )K(y,y 2 )ct(x,dy)('V x h(x,£)(dx))c t {Z,dy 2 )d£, 

1 JxxyJxxy 

- 77 / I g(x,y)K(y,y 2 )ct(x,dy)(V x [h(x,£) + h(£,x)\ (dx)) q(£, dy 2 )d£. (36) 

1 JxxyJxxy 


This characterisation is however limited to functions with compact support in the interior of X. It 
can only give information about how a solution changes within X, it says nothing about what might 
happen on Ti n . Including the boundary terms in the integration by parts/Gauss Theorem used for 
Proposition [53] yields the following additional terms. That these are the correct additional terms is 
part of the assertion of Proposition [57] 
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Definition 54. Let c £ L°° ([0, oo), A4o,oo) and define a vector measure Jfic} on X x y by 


= ( g,Jt[c]) + 


Jr in xy 

1 

2 


g(x,y) ■ n(x ) 


d Ibdryjt, x, dy) 
dt ut(x) ■ n(x) 


+ I int {t,x,dy) - V r • 


u t (x)Ibdvy(t,x,dy) 

ut(x ) • n(x) 


dx 


/ [ K(y,V 2 )g{x,y) ■ n{x) ((h(x,x 2 ) + h{x 2 ,x))) Ihd J^ d '^ Ct (x 2 , dy 2 )dx 2 dx 

l*xyJr in xy u t {x) ■ n{x) 


+ 


'r in xy 


g(x, y) ■ V r / b dry(t, dy)dx - 


f -^bdry (i, 37, dy) 

lr in xy M x ) ■ n ( x ) 


g{x,y) ■ {V r n(x)) T u t {x)dx 


for g £ F 0 * and where Vp is t/ie derivative restricted to directions perpendicular to n(x). Under the 
assumptions on X set out above Vr = (0, ..., 

Definition 55. For c £ L°° ([0, oo), A4o,oo) define time dependent linear operators Gfic] on A4 d by 
(g, G t [c]v^ = -{g- Vu t , i/) + ^t[c]y, /or all g £ B d . 

One can now compactly rewrite the equation from Proposition 1531 as (compare ([3]) ) 


^ ( g , d) = (fiS ^t) + (g, + ( 9 , Jt[c]) (37) 

for all g £ B d such that V • g £ D. The additional terms introduced in Definition [55] are not seen by 
the smaller class of test functions used in Proposition [531 


Proposition 56. Assume HI holds and that c € L°° ([0, 00 ), Ado, 00 ), then there is a strongly contin¬ 
uous, bounded propagator P* ,s [c] on A4 d with 


< ed^coCfillclll^+IIIVulU) min(t-s,t 0 ) 
M d ->M d ~ 

and for f £ D, g, £ M d 



L (/, V‘-’{c]fi = (/, (U t + G,[c]) . 


Proof. This follows the same perturbation argument as Proposition [23] since by duality 

Gt[c] 


M d ^M d 


<|if 0O C' 1 |||c||| fl * + |||V«|||. 


□ 


Proposition 57. Let 13 hold; assume further that either d = 1, H2 holds and inf t,x u t{x) > 0 or H3 
holds for general d; assume further that c £ L°° ([0, 00 ), Mo, 00 ), then (GZD has a unique solution 


v t = V 


t, 0 


[c\v 0 + J P 4,s [c]J s [c]ds € C ([0,oo), ( y M d , ||-|| B *)) 


with initial condition vq. This solution is in L°° ([0, 00 ), 00 ) provided uq £ Adg^ and thus (iden¬ 

tifying the measure with its X-density) also in L°° ([0, 00 ) x X, At(T)xv) • 

Proof. Existence and uniqueness are immediate for this linear problem. Continuity in the F*-norm 
follows from the strong continuity in t of V t,s . 

That the propagators V t,s preserve A4g ^ can be seen by analogy with Proposition[5TJ DefinitionloTl 
expresses the Jt[c] as a sum of Jt[c] and a term concentrated on the inflow boundary. Under H3 
Jt[c ] £ M.q oo and so one argues as in Propositions IT3M441 to show that jjj U t,s [c] J s [c]ds has a density 
with respect to Lebesgue measure on X. 

In the case when only H2 holds, then the x-derivatives of h in (13611 may only exist in a distributional 
sense. However, under H2, the measure V x /i(x,£)( dx) can be expressed as a sum of an absolutely 
continuous part with a bounded density and a finite number of atoms a^(t)S clk with a^ft) £ K, 
ak € X C R. When d = 1 each of these atoms is like a simpler version of the boundary part of the 
inception measure, which in this case reduces under the assumption 12 (see §E2>D to 7 b dr y (C dy)do(dx). 
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The boundedness (uniform in t and k) of the ctk(t) is immediate from the boundedness of K and c 
and since u is bounded away from 0 the analysis of Propositions I73R LH1 applies to show that for each 
k and all t 

[ P i ’ s [c]a fc (t)<5 afc ds € Mq :00 (38) 

Jo 

with a global in time bound in the Ado,oo- n orm. □ 

Proof of Theorem [73]. For d = 1 assume without loss of generality that X = [0, L ) for some L > 0 and 
Tin = {0}. The boundary condition c(t,0, dy) = Jbdr ^*Q°) d ^ is given by Proposition [50] The presumed 
derivative v from Proposition [57] is then used to construct 

c(t,x,dy) = Jbdr y^’ 0 ’ dj/ ) + f v(t,£,dy)d£, (39) 

M o) Jo 

which is readily seen to be a strong solution to ([2]) and therefore to be in the same L°° ([0, oo), A4o,oo) 
equivalence class as c. Therefore (a version of) c is in L°° ([0, oo), VF 1 ’ 00 (X, jM(3^)tv)) and since 4rc 
can be expressed in terms of c and the result follows. 

This argument does not generalise easily to more then one space dimension. However the existence 
of a weak derivative was shown in Proposition 01J] and under H3 Proposition [37] shows that this weak 
derivative in fact has an L°° density. The boundary condition comes from Proposition [50] □ 


5 Discussion 


This paper proves the well posedness of an equation for measures, modelling the creation and coagu¬ 
lation of particles in a flow, for example a flame, for which stochastic approximations were studied in 
14]. In that work the existence of one or more non-negative solutions was proved under somewhat less 
general assumptions there by constructing the solutions as limits of stochastic approximations. The 
present work extends this result by showing that there is in fact only one solution to the equation for 
a given initial condition and thus that all limit points of the approximating sequence from [14] are the 
same and those approximations converge rather than merely having convergent sub-sequences. The 
present work incidentally provides an additional, less constructive proof of the existence of a solution 
to ©. 


It is proved here and in 14] that solutions to ([2]) have a density with respect to Lebesgue measure 
on X and that this is uniformly bounded in time and in X. The differentiability of the density 
is established here even for delocalisations that are of bounded variation, but only in one spatial 
dimension. This result does not extend in full generality to higher spatial dimensions- it is easy to 
imagine two parallel streams of particles that never mix and therefore not even continuity over the 
dividing line in the flow, much less differentiability, is to be expected. It seems therefore likely that 
the discontinuous, cell based delocalisation of the coagulation interaction used for numerical purposes 


m 


15] is not well suited to more than one spatial dimension and that smoother delocalisations should 


be used. Similar methods have been used for the simulation of Boltzmann gases 
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A The flow field 


Definition 58. Let s, t € R and define the flows by 

d 

—® 8jt (x) = u t ($ s ,t(a:)), ®s,si x ) = x - 

$ is a vector, so in more than one dimension it is necessary to distinguish between the matrix 
V<J>, which is the subject of the next two propositions and the divergence, a real number V • <h, which 
occurs in connection with the velocity field u. 


Proposition 59. 

d 

—d> Sit (x) = -V$ s , t {x)u s (x). 

Proof. 


$s,t(x) - x ) 

hm —--- 

<5\o 6 


lim —--- 

<5\o 5 

V$ s flx)(x-<S> s _sA x )) 

= lim--- 

5\o 5 


-V$ s flx)u s (x). (40) 


The right sided limit is dealt with similarly. 


□ 


Proposition 60. 


3 -ll|V«|||(i-«) < ||v$ a)t (x)|| I 


< e l||V«|||(t-«)_ 


and 

det V$ s , t (x) = 

Proof. For the first statement one has 
both smooth, 


det V$t,s(x) = e -fl V - u r(*t,r(x))dr 

ut ($ s ,t(x)) so that, since u and therefore <h are 


d 


V<F S Ax) = V'U) 




(41) 


and the result follows by an application of Gronwall’s inequality. 

The result for the determinant is known as Liouville’s formula. One checks by row operations that 
det Vut (<3? Si t(x)) V<h Si t(x) = Tr(Vrt) det (V<3? Sj t(x)) and the result that follows by solving the resulting 
ODEs. □ 
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